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Abstract 

In the first part of this work we calculate the high frequency magnetoelectric susceptibility of a 
simultaneously ferroelectric and canted antiferromagnetic (also know as weak ferromagnetic) thin 
film with magnetostrictive magnetoelectric coupling. We show that a dynamic coupling exists be- 
tween the ferroelectric and optic antiferromagnetic excitations. In the second part of the paper, 
we calculate using an effective medium method the susceptibility of a heterostructure comprising 
alternating thin films of such a material together with a ferromagnet. Dipolar magnetic fields serve 
to couple the ferromagnetic and optic antiferromagnetic modes, which in turn couples the ferro- 
magnetic and ferroelectric excitations. This provides a mechanism for creating "electromagnon" 
modes in the microwave regime which may be useful for applications. 
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I. INTRODUCTION 



In the 1970s, Bar'yakhtar and Chupis calculated the high frequency magnetic, electric 
and magnetoelectric susceptibility of a model ferroelectric ferromagnet using second quan- 
tization of the electric polarization P and the magnetization M fields They noted 
that the equilibrium directions of P and M must not be parallel or perpendicular in order 
for there to be a dynamic magnetoelectric coupling and in order for the existence of cou- 
pled excitations, known as "electromagnons." Maugin found a similar result [2j. However, 
most simultaneously magnetic and ferroelectric materials (known as multiferroic) are not 
ferromagnetic and have a more complicated spin structure. 

Later, in 1982, Tilley and Scott used a Landau- Ginzburg free energy and equations of 
motion to calculate the full high frequency susceptibility of the antiferromagnetic dielectric 
BaMnF 4 [3j . They were able to explain the observed frequency dependent dielectric anomaly 
in BaMnF 4 by including a magnetoelectric coupling term of the form (/3ip + (32p 2 )M x L z , 
where p is the dielectric polarization, M = M a + M b and L = M a — M b - The subscripts 
a and b denote the two antiferromagnetic sublattices. fti and j3 2 give the strengths of the 
magnetoelectric coupling. This magnetoelectric coupling term is a Dzyaloshinskii-Moriya 
type term {4], Q that causes a canting of the antiferromagnetic sublattices (hence the material 
is called a weak ferromagnet) that may be altered by application of an electric field 6|. 

A similar term has been used to model weak ferromagnetic BiFeC>3, the only known room 
temperature magnetoelectric multiferroic material. deSousa and Moore demonstrated how 
a coupling, P ■ M a x M b , could lead to electric field control of magnon dispersion with 
potential applications to spin wave logic devices Q]. 

This type of coupling in BiFeO"3 now seems unlikely since there is no observation of a 
change in the weak ferromagnetic moment when P is reversed or when an electric field 
is applied s[ However, a magnetoelectric coupling energy which is always symmetry 
allowed is (lOj 

S= (J + TP 2 )M a -M b , (1) 

where J is the antiferromagnetic exchange constant which is perturbed slightly by the in- 
fluence of the electric polarization P. Ionic distortions in the displacive ferroelectric cause 
changes to the effective exchange interaction between spins, giving rise to this magnetostric- 
tive magnetoelectric coupling. We will show that this type of coupling, together with a 
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weak canting of antiferromagnetic sublattices (so that M a /b an d P are not exactly perpen- 
dicular), allows for a dynamic magnetoelectric coupling. It is the "optic" antiferromagnetic 
mode, where the two antiferromagnetic sublattices oscillate out-of-phase, which hybridizes 
with the dielectric mode. 

In Sec. |TT] we calculate the full frequency- dependent magnetoelectric susceptibility tensor 
analytically from a starting free energy for a thin film ferroelectric weak ferromagnet. Com- 
ponents of the susceptibility tensor have poles at the magnetic, electric and magnetoelectric 
resonant frequencies. 

In Sec. II III we extend our calculation to consider a heterostructure containing alternat- 
ing ferromagnetic and ferroelectric weak ferromagnetic layers. Using a particular effective 



medium method for long wavelength dielectric excitations [11] and long wavelength mag- 
netic excitations the susceptibility is found analytically and reduces to known limits. 
We find that the ferromagnetic resonance couples to the optic antiferromagnetic mode via 
dipolar fields and hence also couples to the dielectric mode in the weak ferromagnet. There- 
fore there is a magnetoelectric resonance in the low GHz (microwave) regime, whereas in 
single-phase magnetoelectric materials the resonances are usually all in the infrared regime. 
By combining such a ferromagnet and weak ferromagnet it is possible to create an effective 
material for tuning magnetoelectric response frequencies. Applications exist to microwave 
signal processing using applied electric fields jjj], Q or even to designing left-handed mate- 
rials in small frequency ranges 16] . 



II. WEAK FERROMAGNET SUSCEPTIBILITY 



A. Geometry and energy density 

The geometry of the ferroelectric weak ferromagnet is shown in Fig. [1] The electric 
polarization P lies along the x direction. The two antiferromagnetic sublattices M a /b he 
perpendicular to the electric polarization, predominantly along the y direction, but are 
canted in the z direction by an angle 9. This angle is exaggerated in Fig. [I] and in BiFeOs, 
for example, is given by 9 ~ 0.14° [171 ] . When considering the thin film geometry, the film 
thickness is in the z direction. This minimizes depolarizing plus demagnetizing energies. 
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FIG. 1: The ferroelectric weak ferromagnet geometry. 



The energy density of the antiferromagnet is given by: 

Safm = ( J + T(P x f)M a ■ M b — K [(Ml) 2 + (M y b ) j 

+D(MlM b z - M z a Ml) - h M + 2n(M* + M 6 2 ) 2 . (2) 

The first term represents the antiferro magnetic exchange interaction with J > 0. A weak 
contribution to the exchange constant is due to so-called "isotropic" or magnetostrictive 
magnetoelectric coupling with strength given by T [lOj. It arises since the soft-phonon mode 



associated with the electric polarization in the film P is coupled to the magnetic system 
through magnetostriction. The second term is a uniaxial anisotropy energy which favors the 
alignment of the sublattice magnetization in the y direction. The third term in Eq. (j2J) is 
the Dzyaloshinskii-Moriya interaction with a Dzyaloshinskii-Moriya vector D = Dx giving 
canting in the z direction. The fourth term describes the interaction with a small driving 
field h. Finally, the last term in Eq. (j2J) is the demagnetizing term in CGS form, assuming 
that the thin film geometry has an interface containing both sublattices. 
The energy density of the dielectric part of the system is given by: 

S FE = -l£(p*) 2 + lA(P*) 4 -e.P 

+\u[{P y ? + {P z ) 2 ]+^{P z ?. (3) 

£ and A are phenomenological Landau coefficients giving a spontaneous polarization in the 
x direction. The isotropic magnetoelectric coupling constant V (see Eq. (J2J)) alters £ by 
a small amount. A one-dimensional model for the spontaneous polarization is valid when 
examining small amplitude dynamics about equilibrium. The third term in Eq. ([3$ is the 
interaction of the dielectric with a small driving field e. The fourth term describes the 
strength of the dielectric stiffness £j_ > of the material in the y and z directions. We make 



a simplifying assumption that the system is isotropic in the y-z plane. The last term is the 
depolarizing energy density. 



B. Equations of motion 

From the free energy, the equations of motion for the magnetization and polarization 
can be found using the Landau-Lifshitz (or torque) equation and the Landau-Khalatnikov 
relaxation equation respectively 

dM ljr / d£\ 

^ M >< -Ui7 > ( 4 ) 



dt V dM 

d 2 P „ / dS 



dt* f \dP)' (5) 

where the derivatives, — j^ and — Jp, represent the effective magnetic field and the effective 
electric field acting on the systems. 7 is the gyromagnetic ratio and / is an effective inverse 
mass term for the dielectric oscillations. We ignore damping in both equations. 

The equations of motion are obtained by substituting Eqs. ([2j) and (j3J) into Eqs. (j4j) and 
(jSJ), assuming oscillating solutions that vary in time according to e~ luJt , and then linearizing 
the resulting equations. The linearization is done by splitting the two sublattice magneti- 
zations and the polarization into static and small dynamic parts and then ignoring terms 
which are quadratic in small dynamic terms. If dynamic parts are denoted by lower case 
letters, then according to the geometry shown in Fig. [1] the equations are linearized using: 

M a = (a x , M cos 6 + a^Mo sin 6 + a z ), (6) 
M b = (b x ,-M cos6 + b\M sm6 + b z ), (7) 

p = (p + P x , p y, P z ). (8) 

The equilibrium canting angle 9 is given by minimizing Eq. (j2J): 

6 = — arctan ( = ] , (9) 

and the equilibrium polarization in the x direction, Pq, is given by minimizing Eq. (jHJ) plus 
Eq. ©: 

/(£-2rM n 2 [- cos 2 6 + sin 2 01) , 

P o = V- qL -^ -• (10) 
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a y = a" 



Combining Eqs. ©-(jHD we obtain magnetization equations: 

-a y ([2H d + H ex (P ) + H a ] sin 9 + H DM cos 9) 
-a z ([H d + H ex (P ) + H a ) cos 9 - H DM sin 9) 

+b y (H DM cos 9 + P e:c (P ) sin 6) - b z ([H d + H ex (P )] cos 9 - H DM sin 6) 
-p x 2TP M$ sin(2#) + M h z cos9 - M h y sm6 (11) 
i x ([2H d + H ex (P )} sin 9 + P OM cos 0) - b x H ex (P ) sin 
+M o /i a, sin0 (12) 
' ([H ex (P ) + H a ) cos 9 - H DM sin 0) + b x H ex (P ) cos 
-M h x cos 9 (13) 
-^6* = a y (P ex .(P ) sin 9 + H DM cos 9) + a z ( [P d + P ei .(P )] cos - H DM sin 0) 
-ft 4 ' ([2P d + P e ,(P ) + H a ] sin 9 + H DM cos 0) 
+b z ([H d + P ex (P ) + Pa] cos 6 - H DM sin 9) 

+p x 2FP M$ sin(2#) - M /f cos 9 - M h y sin 9 (14) 
V = -a :c P e:E (Po)sin^ + 6 3; ([2P d + P ex (P )]sin^ + P DM cos^) 

+M h x sin (15) 



a = a ' 

7 

iuj . 
1 
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- = -a x H ex cos 9 - b x ( (P ) + P a ] cos - H DM sin 5) 

+M o /i a: cos0, (16) 

where the effective exchange, anisotropy, demagnetizing and Dzyaloshinskii-Moriya magnetic 
fields are given respectively by H ex (P ) = M (J + TP 2 ), H a = 2KM , H d = 4vrM and 
Hdm = MqD. We write H £X (Pq) as H ex below to shorten the notation. 

The only component of the dielectric polarization to couple with the magnetization equa- 
tions is p x . It's equation of motion is given by: 

- —f = (£ - 2TMl\- cos 2 9 + sin 2 9] - 3AP 2 ) p x + e x 

-2TP M (-a y cos 9 + b y cos 9 + a z sin 9 + b z sin 9) . (17) 

It can be seen that if the canting were to vanish, then 9, a y , b y — > 0, and the magnetic and di- 
electric equations of motion would not be coupled. So although the magnetoelectric coupling 
enters into the exchange interaction, rather than the Dzyaloshinskii-Moriya interaction, it 
results in a dynamic magnetoelectric coupling. 



The equation of motion for the remaining two components of the dielectric polarization 
are: 

- jp y = -UP v + e y (18) 
-jP z = -{U + ^)p z + e z . (19) 

C. Susceptibility 

The seven equations of motion Eqs. ffTTl) - (fTTj) can be used to solve for 
{a x ,a y ,a z ,b x ,b v ,b z ,p x } analytically as a function of driving fields h and e x . First we set 
h z and h x = h y = e x = 0. Then the equations for a and b are symmetric under the 
transformation b x —>■ —a x , b y — > —a y and b z — > a z . This is the so-called "optic" antiferromag- 
netic mode where the two antiferromagnetic sublattices oscillate out-of-phase. Eqs. (TTTT)- HTTT) 
reduce to: 

%UJ -a x = -a y ([2H d + 2H ex + H a ] sin 9 + 2H DM cos 9) + M h z cos 9 



7 



-a z ([2H d + 2H ex + H a ] cos 9 - 2H DM sin 9) - p x 2TP M£ sin(2#) (20) 



a y = a x {[2H d + 2H ex \sin9 + H DM cos6) (21) 

7 

-—a z = a x (H a cos 6 - H DM sin 6) (22) 

7 

p x = - 2TMl [- cos 2 9 + sin 2 9] - 3 AP 2 ) p x 

-4TP M (-a y cos9 + a z sm9). (23) 



L0 2 



f 



Eqs. (!2T!)- (!23l can be substituted into Eq. (|20|) to get an equation involving only a x : 

a x = ^ M « hZ ™\ (24 ) 

CO 2 — LO Z 

where the optic antiferromagnetic frequency u Q is given by: 

( [2H d + 2H ex ] sin 9 + H DM cos 9) ( [2H d + 2H ex + H a ] sin 9 + 2H DM cos 9) 



w2 



7^ 

+ (H a cos - H dm sin 0) ([2i?d + 2i7 e:E + H a ] cos - 2H DM sin i 
_ 8r 2 P 2 M 3 sin(2#) (cos 9 sin [2H d + 2H ex + H a ] - H DM ) 
2rM 2 [- cos 2 9 + sin 2 9} - 3AP 2 ) ' 



(25) 
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Since the Dzyaloshinskii-Moriya canting angle is small, u is approximated very accurately 
by taking the limit sin 9 — > and cos 9 — > 1 . This gives: 

4 ~ 2H 2 DM + H a (2H d + 2H ex + H a ) . (26) 

Ignoring the effective Dzyaloshinskii-Moriya field H^m, this f req uency agrees with the well- 
known result for thin film antiferromagnets with no canting [18j |. Eq. (!26l) also agrees with 
the resonant frequency calculated previously for bulk canted antiferromagnets when instead 



H d = 



2Q| 



The xz component of the magnetic susceptibility x™ z = {a x +b x ) / h z is zero since a x = —b x . 
Similarly, is zero. The nonzero susceptibility components due to h z are xTz = {a z + b z )/h z 
and the electromagnetic susceptibility xt™ = p x /h z which are given exactly by: 

-2 7 2 M cos 9 {H a cos 9 - H DM sin 9) 
Xzz oo 2 -uj 2 ' 1 j 

em 4/rP o7 2 M 2 cos 9 (cos 9 sin 9 [2H d + 2H ex - H a ] - H DM ) 
XxZ ■ W 

X x ™ has a pole at the optic antiferromagnetic mode frequency u Q and also at the ferroelectric 
mode frequency Uf e given by 

^ = -f + 2rM 2 [- cos 2 9 + sin 2 9} + 3AP 2 . (29) 

Driving fields e x excite the same modes: the ferroelectric mode and the optic mode. 
Eqs. ( I20l) -(j23l) remain the same apart from the removal of h z and the inclusion of e x . Fol- 
lowing the same working, it is found that xTx = ( flZ + b z )/e x = xTzi which is given in 
Eq. (|28]) . The only other nonzero component appearing due to e x is: 



-/{^ 2 - 7 2 {[2H d + 2H ex ] sin 9 + H DM cos 9) {[2H d + 2H ex + H a ] sin 9 + 2H DM cos 9) 
- 7 2 (H a cos 9 - H DM sin 9) ([2H d + 2H ex + H a ] cos 9 - 2H DM sin 9) } 

? e )(- 2 -- 2 ) (30) 




2 2 \ / 2 



^ 2 - ^ 2 e ' 



Next we set /i x ' 7^ and h z = h y = e x = in Eqs. (fTTl) - (IT7j) to solve for the susceptibility 
components x%x (i — x,y,z). The equations for a and b are symmetric under transform of 
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a x — > b x , a y — > b y and a z — » — 5 Z , which corresponds to the antiferromagnetic sublattices 
oscillating in-phase and is referred to as the "acoustic" mode. The magnetoelectric coupling 
term at the end of Eq. (1171) vanishes and hence p x = when the system is driven by a 
magnetic field in the x direction. 
Eqs. dHD-(HZD reduce to: 

- — a x = -a y ([2H d + H a ] sin9) - a z H a cos 9 (31) 

7 

-— a y = a x (2H d sin 9 + H DM cos 0) + M sin fl/i* (32) 

7 

-^a z = a x ([2^ + H a ] cos - H DM sin 0) - M cos 9h x . (33) 

7 

The two nonzero components of the susceptibility xTx an d X™ x are given by: 
m 2 7 2 M (sin 2 9 (2H d + H a ) - cos 2 9H a ) 



a 



v 



i2 7 Mo (u; 2 sin g - 7 2 ff a cos g ( [2H d + 2if e:E + H a ] cos g sin g - H DM [1-2 cos 2 g] )) 
" yx = (cu 2 - ul) ^ 

where the acoustic antiferromagnetic resonant frequency is 

4 = H a (2H ex + H a ) cos 2 9 + # d # DM sin(2#) + 2# d (2# d + H a ) sin 2 0. (36) 

7 2 

Once again, if we make the approximation sin# — ► 0, then this expression reduces to the 

nnri 

known acoustic frequency (whether there is canting or not) given by [181 . 1 191 . [20| 



^~H a (2H ex + H a ). (37) 

Making h y ^ and h x = h z = e x = 0, we find that driving fields in the y direction also 
excite the acoustic antiferromagnet mode and do not excite a dielectric mode. By symmetry, 
the susceptibility component Xx]j = ~X y n x an d so has already been found. The remaining 
susceptibility component \™ y is found to be: 

m 2-f 2 M sm9(2H d sm9 + H DM cos9) 

Xyy - 2 2 • ^°>> 

This component vanishes as 9 — > since then the linearization of the antiferromagnetic 
sublattice magnetizations requires that there is no dynamic magnetization in the y direction. 

Finally, there are two non-zero components of the electric susceptibility given by exam- 
ining Eqs. (fig]) - (113]) : 

f_ 

-uj 2 + 



Xty = — rWT- (39) 



x * z = -^ 2 + /(a + 47r)- (40) 



The total susceptibility tensor for the ferroelectric weak ferromagnet geometry takes the 
form: 



(mm q q 



X 



^xx /Vxy 
A-yx X-yy 







q q m v me 






V 



lzi 










\ 











(41) 



" Axz Am 
X e yy 

xlj 

with the eight independent components being given in Eqs. (1271) . fl28l) . fl3Tl) . (1341) . fl35l) . (1381) . 
(|39]) and gQJ. 



III. WEAK FERROMAGNET/FERROMAGNET HETEROSTRUCTURE 
A. Energy density and equations of motion 

We now consider a heterostructure comprised of alternating thin films of a ferroelectric 
weak ferromagnet, as illustrated in Fig.fU with thickness d w and a ferromagnet with thickness 
df in the ^-direction. We shall assume that the ferromagnet has a uniaxial anisotropy in 
the y direction with strength Kf and dielectric stiffness components given by a x , a y and a z . 
Then the energy density of the ferromagnetic film is given by: 

g FM = -K f (Mj) 2 -H r M f + l - J2 «*b}) 2 - E f ■ Pfl (42) 

i=x,y,z 

where Mf = (m x ,Mf + m y ,m z ) is the linearized magnetization and is the dynamic 
dielectric polarization and so is denoted by a lower case letter. H f and Ef are the dipolar 
magnetic and electric fields respectively. They have been written in upper case to emphasize 
that these may have a static as well as a dynamic part. 

We rewrite the energy density for the ferroelectric weak ferromagnet shown in Eqs. ([2]) 
and ([3]) so that the thin film demagnetizing and depolarizing terms are discarded and the 
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dipolar fields are written in a corresponding way as to in the ferromagnet: 

<&afm = (J + T(P x ) 2 )M a ■ M b — K [(My) 2 + (M y h f 

+D(M*M§ - M z a M'l) — H w M (43) 
See = -\t {P x f + \ A (P x )' ~ E w P 

+~t± [{P y f + (P z ) 2 } • (44) 

To calculate the analytic susceptibility and the resonant k = frequencies of the het- 
erostructure, we use an effective medium method which requires that Maxwell's boundar 
conditions for dipole fields are satisfied at the interfaces between the materials [ll 
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This method can also be used to numerically calculate the frequencies of long wavelength spin 
waves with finite wavevector (k ^ 0) in an approach known as entire-cell effective medium 
method 2l|, |22| and gives results in good agreement with more computational-demanding 
methods for including dipolar interactions. 

Maxwell's boundary conditions relate the dipolar fields in the ferromagnet (Hf, Bf — 
H J — AnMf, Ef and Df = Ef + 4nPf) and the weak ferromagnet (H w , B w = H w + 
47r(M a + M b ), E w and D w = E w + 4ttP w ) according to: 

H* = H x w = h x (45) 

H) = H* = h" (46) 

Hf + AnM z f = H z w + Att(M z + M z h ) = C (47) 

E x = E x w = e x (48) 

E) = El = e y (49) 

E z f + AixP z = E z w + 4nP z = D, (50) 

for the geometry shown in Fig. [TJ The constants C and D are defined for ease of notation 
in what follows. In particular, the out-of-plane (z) boundary conditions couple the dipolar 
fields to the magnetization and electric polarization in both materials. For this reason we 
must calculate the iz components of the susceptibility tensor first in order to properly take 
into account dipolar effects. 

All of the dipolar fields are in fact dynamic, apart from since from the linearization 
Eqs. and ([7]) together with the boundary condition Eq. fj4T|) 

H z w = C- 8rrM sin 9 - An(a z + V). (51) 
11 



It is the dynamic part of if*, namely h z w = C — 47r(a z + b z ), which drives the magnetization 
and which we need to find in order to calculate the dynamic effective medium susceptibility. 

Substituting the energy densities Eqs. and boundary conditions Eqs. (r4"5l) - (l5T)j) 

into the equations of motion Eqs. (jlj) and §5§ we obtain the following equations. For the 
ferromagnet we have: 

-^=-m^, + H, ) + M f C (52) 
- % ^-m z = m x H af + M f h x , (53) 

where the effective anisotropy and static dipolar fields in the ferromagnet are given by 
H a f = 2KfMf and = 4nMf respectively. We assume that the gyromagnetic ratio 7 is 
the same for both the ferromagnet and the weak ferromagnet. 

The equations of motion for a and b are the same as shown for the weak ferromagnet 
thin film in Eqs. (|TT]) - (fl6|) apart from the replacement of h z — > C. This is deceptive as it 
appears that the boundary conditions have simply created a thin film demagnetizing effect. 
However, since the driving field is h z w = C — 4ir(a z + b z ) rather than C, this is not the case 
as will be shown later. 

The linearized electric equations of motion are: 

- jV x f = -a x p x f + e x (54) 
~Y pV f = ~ a y pV f + eV ( 55 ) 

-jV) = -{a z + A7r)p z f + D (56) 
- J-Pl = (e - 2TM 2 [- cos 2 6 + sin 2 9} - 2K ME M 2 - 3AP 2 ) f w + e x 



-2TP M 2 {-a y cos 9 + b y cos 9 + a z sin 9 + b z sin 9) , (57) 



-rVl = -UpI + e y (58) 

Ju 



I w 

,2 



-7-K = ~(C± + ^)P Z W + D, (59) 

Jw 



where ff and /„, are the effective inverse mass terms of the respective dielectric materials. 
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B. Effective medium susceptiility 



As shown in Sec. [TTJ the magnetizations and electric polarizations can be found as a 
function of dipolar field and then the susceptibility can be derived. As already mentioned, the 
Xiz (i = x, y, z) components of the susceptibility must be found first for the heterostructure 
which involves setting C ^ and D / and ignoring all other dipolar field components. 
In the effective medium approximation the fields in the two materials are averaged (for 
example, (m x ) = dfm x + d w (a x + b x )) giving susceptibility components: 

m = d f m ' 1 + d ™( al + bi ) = W (m) 

Xiz d f {C - ATim z ) + d w (C-Ana z - 4ir¥) (h z ) ' 1 ' 

d f p) + d w pj _ (p») 

Xiz d f (D - Anpj) + d w (D- Anpzj ~ (e*> ' 1 ) 

Y em = d fP} + d w Pj = M ( 62 ) 

Xtz d f (C-Anm z ) + d w (C-4na z - 4nb z ) (/i 2 )' 1 ' 

Xiz d f (D - Anpj) + d w (D - 4npz w ) ~ (e 2 ) ' 1 ) 

Each term is weighted by the corresponding film thickness d w or df. 

The other susceptibility components can then be found. For example, by setting h x ^ 
the Xii an d XiT components can be found according to: 

d f m* + d w (a* + V) - xZ(h z ) - X™ e (e z ) 



Y m 
A.ix 



(d f + djh x 

(™ l )-X? z (h z )-xZ e (e z ) 
(h x ) 



(64) 



em = (f)-x!T(h z )- X Ue z ) 

Aix (h x ) ' 



where the susceptibility components on the right hand side of Eqs. ( 1641) and ( 1651) were found 
in the previous step using Eqs. (|60l)-(l63l). 

It should be noted that the method described is identical to existing effective medium 



methods 
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21 



22|, apart from the fact that we treat both dielectric and magnetic 
systems for the first time. This only works since the equation of motion for p z w and p z ^ 
are not coupled to the equation of motion for m z , a z and b z . If the out-of-plane dielectric 
and magnetic oscillations were coupled, then the effective dipolar fields (h z ) and (e z ) would 
be coupled and much more complicated expressions for the susceptibility components, as 
compared with Eqs. (160]) - (163]) . would need to be found. This represents a new extension to 
the effective medium method which will be discussed in a later paper. 
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Without providing working, the result for the non-zero components of the frequency- 
dependent susceptibility is 

-d n 2 M f H af (co 2 - cop - d wl 2 M cos 6 (H a cos 9 - H DM sin 6) {co 2 - co 2 ) 
Xzz ~ d f (co 2 - co 2 Q ) {co 2 - co%) + d w (co 2 - col) " wj) 

-i^M f d f {co 2 -cop 

XZ ' " d f (co 2 - cop {to 2 - col) + d ™ V - O - "/) 

i2jM d w (co 2 sin 9 - -? 2 H a cos 9 ([2H d + 2H ex + H a ] cos 9 sin 9 - H DM \l - 2 cos 2 AM 



1 U->: 



co(d f + dj (co 2 - col) 
d f M fl 2 {d f H af (co 2 - cop + d w (H af + 4nM f ) (co 2 - coj)} 
(d f + d w ) {d f (co 2 - cop {co 2 - co%) + d w (co 2 - col) (^ 2 - w?) } 

2d wl 2 (M sin 2 9(2H d + H a ) - M H a cos 2 9) 
+ (d f + d w ) (co 2 - coP 

2d w -f 2 M sin 9 (2H d sin 9 + H DM COS I 
(d f + d w ) (co 2 - col 



(69) 

m _ *u w y imp sun/ ^n d aiu.u ~r upM^uj 



em -d^T^^PpM 2 cos 9 (cos 9 sin9[2H d + 2H ex - H a ] + H DM ) {co 2 - co 2 ) 

XxZ " {u; 2 -co 2 e ){d f (u 2 -coP{co 2 -u%)+d w (co 2 -col){u 2 -co 2 )} ( j 

em tcod f d w ATT 7 HTf w P M 2 M f cos 9 (cos 9 sin fl[2ff d + 2H ex - flj + ff JA /) 

X ^ (d, + O {co 2 - co%) {d f (co 2 - col) (" 2 - «JJ + (c 2 - a&) (c 2 - ^ 2 ) } 1 J 

<*/ {-u 2 /fw + U + 4?r) + rf w {-u 2 /f f + a z + 4tt) 

V d f (co 2 /f f - a z ) (co 2 /f w -U- 4tt) + d w (co 2 /f w - C±) (co 2 /f f - a z - An) 

XL = 7^-Vr(^4- + ^H (74) 



(d f + d w ) \co 2 - f f a a 



co'- — CO 



f e . 



x w (a _lw ^ (,.,2_f.„ + ,.,2 x ^ ) ' ( 75 ) 



-i f <*/// , gj, 

(dy + dj \C0 2 - ffOty ^ - fintl-. 

where the optic antiferromagnetic frequency in a weak ferromagnetic thin film co is given in 
Eq. (125p . the acoustic antiferromagnetic frequency co a is given in Eq. ( 1361) and the ferroelectric 
mode frequency u;j e is given in Eq. (|2"9~|) . In addition, we have new frequencies for the 
ferromagnet in thin film coj and in bulk co^\ 



^ = "a, (Ha, + i*M t ) , (76) 

^ = Hi,. (77) 

An additional frequency associated with the optic antiferromagnetic mode in bulk is given 
by 

2 2 

^ = ^-2H d cos 9 (H a cos A - H DM sin A) . (78) 
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Compared with the weak ferromagnet in isolation (Sec. [TT]) two extra components are 
non-zero in the susceptibility tensor, namely xS, and Xxx- X^z appears since it is non-zero 
in the ferromagnet. XxH arises purely due to the coupling of m x in the ferromagnet to the 
out-of-plane dipolar field and is given by: 



cm 



em _ dfAirm x xz 

(d f + djh*' 1 ] 



C. Limiting cases 



We consider some limiting cases to test the effective medium method. We use xTz ( see 
Eq. dSS}) to demonstrate the results. 

First we consider replacing the ferromagnet with a nonmagnetic material (Mf — > 0). The 
component becomes 

m -d w ^ 2 M Q cos 6 (gg cos 6 - H DM sin 6) 
Xzz d f (^-^ ) + d w (^-u;l) ■ { > 

Then taking the limit that the weak ferromagnetic films are much thinner than the nonmag- 
netic spacers (dj » d w ), the isolated thin film result from Sec. [TT] is recovered (Eq. ( |27j) ). 
namely 

mfiim cos 9 (H a cos 9 Hj->m sin 9) 
Xzz -r—, • (81) 

Next we consider removing the ferromagnetic layers (df — * 0). This gives, 

mbuik -27 2 M cos 9 (H a cos 9 - H DM sin 9) 

= ~2 ~2 ^ 2 > 

which is the same as the result found for the isolated thin film weak ferromagnet (Eq. (1271) ). 
apart from the pole being at the bulk frequency u m rather than at the thin film frequency u . 
Hence the bulk limit is correctly recovered. Similarly, the bulk ferromagnetic susceptibility 
is recovered in the limit of the weak ferromagnet vanishing: 

m.bulk -1 2 H af M f 
Xzz - — 2 "a—- (M) 

The effective medium method recovers the correct limits for both bulk and thin film 
geometries and therefore seems reliable. 
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D. Dynamic magnetoelectric coupling 



What is most significant when examining Eqs. (!66|) - (!75|) is that the effective medium 
susceptibility is not, in general, given by an average of the susceptibility in each film. This 
means that instead of finding poles in xTz a ^ the ferromagnetic bulk frequency oj^ and at 
the optic antiferromagnetic bulk frequency u m) we find two resonant frequencies given by 
the solution to: 

= d f {uo 2 - uol) {uo 2 - uo%) + d w {uo 2 - ool) (uj 2 - u 2 f ) . (84) 

These two frequencies correspond to modes that are common to both materials and are a 
signature of the fact that the out-of-plane dipolar magnetic fields serve to hybridize the 
ferromagnetic and optic antiferromagnetic resonances. In Sec. [Ill we showed that for a 
weak ferromagnet with magnetostrictive magnetoelectric coupling, the ferroelectric and optic 
antiferromagnetic modes are coupled. This in turn means that the ferromagnetic resonance 
is coupled to the ferroelectric mode. Examining the expression for y^™ in Eq. (17TT) . it is 
indeed seen that there is a magnetoelectric resonance involving the ferromagnet. 

To demonstrate that one of the magnetoelectric resonant frequency may be in the GHz 
regime through this indirect coupling of ferroelectric-optic-ferromagnetic modes, approxi- 
mate frequencies for a NiFe/BiFeO"3 (ferromagnet /weak ferromagnet) heterostructure are 
calculated. Equal volumes of both materials are assumed [d^ = d w ). The relevant frequen- 
cies of the isolated films and bulk samples are given in Table [H Substituting these into 
Eq. (I84l) . we find two of the three magnetoelectric resonant frequencies in the heterostruc- 
ture at 4.07 GHz and 548.0 GHz. The former value shows how such a heterostructure may 
be designed to give dynamic magnetoelectric coupling in the microwave regime. With a 
change in the ferromagnet used, application of an applied magnetic field, or a change in the 
relative thicknesses of the two materials, this frequency can be tuned. 

We should stress that the only mechanism in this model for a dynamic magnetoelectric 
coupling between the weak ferromagnet and the ferromagnet is through dipolar fields. In a 
real system exchange coupling at the film interfaces may also lead to a dynamic magnetoelec- 
tric coupling by coupling the ferromagnetic and optic antiferromagnetic modes. Exchange 
coupling leads to an asymmetry between the two antiferromagnetic sublattices and so the 
effective medium susceptibility must be found numerically rather than analytically. Also, 
for relatively thick films, the exchange coupling only represents a small contribution to 
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TABLE I: The resonant frequencies of NiFe (thin film and bulk ferromagnetic modes) and BiFeOs 
(thin film and bulk optic modes). These are estimated by assuming 7 = 2ttx 2.8 x 10 6 Hz/Oe, that 



2$,H ex = 2.7 xlO 5 Oe 



for NiFe H af = 10 Oe and M f = 867 Oe, and that for BiFe0 3 , H a = 880 Oe 

24], Mo = 750 Oe [3] and H DM = 1400 Oe. The value for H DM is estimated using the canting 
angle 6 = 0.14° {3] together with Eq. ©. 







UJ 


U>o* 


5.81 


0.176 


552.6 


543.3 



the energy density and so will not change the resonant frequencies significantly from those 
calculated here. 



IV. CONCLUSION 



In this work we have shown that a magnetostrictive magnetoelectric coupling together 
with a canting of antiferromagnetic sublattices (known as weak ferromagnetism) in a ma- 
terial leads to a dynamic coupling between ferroelectric and optic antiferromagnetic excita- 
tions. Such a model is applicable to known multiferroic materials, such as BiFeOs. Hybrid 
magnetoelectric excitation (or electromagnons) are interesting for probing the origin and 
strength of magnetoelectric coupling but also may have application to high frequency signal 
processing. Most antiferromagnetic and ferroelectric resonant frequencies are in the infrared 
regime. 

In the second part of the work, we used an existing effective medium method to calcu- 
late the high frequency susceptibility of a ferroelectric weak ferromagnet /ferromagnet het- 
erostructure. This is the first time dielectric and magnetic susceptibilities have been found 
simultaneously using this method. The main result is that the magnetic dipolar coupling 
between the films mediates a dynamic coupling between the ferromagnetic and ferroelectric 
modes. This means that there is an electromagnon in the low GHz or microwave regime. 
Heterostructures may be designed to produce electromagnons in a desired frequency range. 

The strength of the dynamic magnetoelectric coupling is in general weak via this mech- 
anism. In fact, for applications it appears that magnetostrictive/piezoelectric composites 
with an interface strain-mediated coupling are much more promising since they have mag- 
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netoelectric coupling strengths up to 100 times larger than in single-phase materials (see, 
for example, the review article by Nan et al. [26!]). Such heterostructures may be similarly 
treated using an effective medium method, with an appropriately chosen magnetoelectric 
coupling between films. A microscopic entire-cell effective medium method may prove more 
useful for calculating the susceptibility since the magnetoelectric coupling is an interface 
effect. 
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